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ABSTRACT 

5" 

In this letter we study the observed distributions of rotational velocity in a 
en sample of more than 16,000 nearby F and G dwarf stars, magnitude complete 

Oh 1 and presenting high precision Vsini measurements. We show that the velocity 

distributions cannot be fitted by a maxwellian. On the other hand, an analysis 
^ 1 based on both Tsallis and Kaniadakis power-law statistics is by far the most 

appropriate statistics and give a very good fit. It is also shown that single and 
binary stars have similar rotational distributions. This is the first time, to our 
£> ■ knowledge, that these two new statistics are tested for the rotation of such a 

large sample of stars, pointing solidly to a solution of the puzzling problem on 
OO ■ the function governing the distribution of stellar rotational velocity 

cn 

O ' Subject headings: stars: rotation-statistics-evolution-fundamental parameters 
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1. Introduction 

Rotation is a fundamental observable for the study of evolution of stars, providing also 
valuable informations on stellar magnetism, mixing of chemical abundance in stellar interior 
and tidal interaction in close binaries. In addition, it can be argued that the value of the 
rotational velocity of a star at a given evolutionary stage may be a measure of its original 
angular momentum. In that case, the distribution of rotational velocity may also be used 
to study some of the characteristics of the physical processes controlling star formation. In 
spite of a large observational effort carried out by different authors along the past 50 years, 
the problem concerning the nature of the statistical law controlling the distribution of stellar 
rotational velocity has turned into a very puzzling question. On the other hand, there is 
widespread recognition that stellar rotation axes have a random orientation (Struve 1945). 
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The first to derive analytically the distribution of stellar projected rotational velocity, on 
the basis of a Gaussian distribution were Chandrasekhar and Munch (1950). These authors 
first assumed a parametric form for a function f(v), where v is the true rotational velocity, 
then computed the corresponding distribution of the projected rotational velocity Vsini and 
finally adjusted a set of stellar parameters to reproduce the Vsini measurements. Although 
Deutsch (1970) has claimed that the distribution of stellar rotational velocities follows a 
Maxwellian-Boltzmann law, a number of studies have shown a clear discrepancy between 
theory and observations, where observed distributions are not well fitted by a Gaussian or 
Maxwellian function (e.g.: Wolff et al. 1982; de Medeiros et al. 1996). More recently, Soares 
et al. (2006) have shown, based on an analysis of the rotation of low-mass stars in the 
Pleiades open cluster, that the question of the nature of the distribution of stellar rotational 
velocity is not simply a question of which mathematical function model is used. It depends 
on the statistical mechanics applied. 

Here, we investigate the effects of powerlaw statistics on the observed distribution of 
projected rotational velocity measurements of a sample of more than 16,000 nearby F and G 
dwarf stars from the Nordstrom et al. (2004) catalogue. We show clearly that the velocity 
distributions cannot be fitted by a maxwellian. Our analysis is based on both Tsallis and 
Kaniadakis non-gaussian statistics and it is shown that these are by far the most appropriate 
statistics, giving a very good fit. We also show that single and binary stars have essentially 
similar rotational distributions. 

This Letter is organized as follows. In Sec II, based on the formalism presented by 
Deutsch (1970), we present a generalization of the rotational velocity distribution in the 
spirit of Tsallis and Kaniadakis statistics. A brief discussion on the stellar sample is made in 
Sec. III. Our main results are discussed in Sec. IV and we summarize the main conclusions 
in Sec.V. 



During the last two decades a great interest has arisen in the non-extensive Tsallis 
statistical mechanics (Gell-Mann & Tsallis 2004) and, more recently, on the Kaniadakis 
extensive generalized power-law statistics (Kaniadakis 2002; 2005) motivated by various 
restrictions to the applicability of classical, extensive statistical mechanics. 

From the mathematical point of view, Tsallis statistics is based on the g-exponential 
and g-logarithm functions, which are defined by 
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Tsallis and Kaniadakis distribution functions 
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whereas the g-entropy associated with the g-statistics is given by (Gell-Mann & Tsallis 2004) 

S * = -J d 3 pf\n q f = -(\n q (f)}. (3) 

The expressions above reduces to the standard results in the limit q — 1. 

Similarly, the K-framework is based on ^-exponential and ^-logarithm functions, defined 

as 

exp K (/) = (v / l + ^ 2 / 2 + ^) 1/K , (4) 

£K r — re 

M/) = ( 5 ) 

The K-entropy is given by (Kaniadakis 2002; 2005) 

S K = ~ J d 3 pf\n K f = -(\n K (f)). (6) 

Again, the standard results are attained in the limit k = 0. 

The Tsallis statistics has been investigated in a wide range of problems in physics 0. 
In the astrophysical domain, the first applications of this powerlaw statistics studied stellar 
polytropes (Plastino & Plastino 1993) and the peculiar velocity function of galaxy clusters 
(Lavagno 1998). More recently, Kaniadakis statistics (Kaniadakis 2002; 2005) has also been 
studied in the theoretical and experimental context, however the first application with a 
possible connection with astrophysical system has been the simulation in relativistic plasmas 
(Lapenta et al. 2008). More recently, by considering the distribution of stellar rotational 
velocity for low-mass stars in the Pleiades open cluster, Carvalho et al. (2008) showed that 
k and q distributions gives a good fit for the observed distribution. 

Following the model given by the authors (Soares et al. 2006 and Carvalho et al. 2008), 
it is possible to show that the Kaniadakis and Tsallis distributions are given by 

F(Q) = exp K (^Pj , F(Q) = ex Pq (^Pj (7) 

where Q is the non-dimensional quantity wj, j being a parameter with the dimension oj^ 1 . 

Here, it is worth mentioning that the standard distribution of the true rotational velocity 
V for a star sample is F(V) ~ V 2 exp(—V 2 ). As shown by Deutsch (1970), the standard 



1 For a complete and updated list of refences sec http://tsallis.cat.cbpf.br/biblio.htm 
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observed distribution of the projected rotational velocity Vsini, for a random orientation 
of axes, must be given by 4>(y) ~ yexp(— y 2 ) (Kraft 1970), with y = Vsini. Henceforth, 
the ^-distributions <j> K (y) (g-distributions 4> q {y)) should reproduce the standard one, in the 
same way as F K (v)(F q (v)) recovers F(V) in the k — (q — 1) limiting case. Therefore, 
by considering these arguments, we introduce the following distribution function for the 
observed stellar rotational velocities 



The rotational velocities Vsini used in the present analysis were taken from the ro- 
tational all-sky, magnitude-limited and kinematically unbiased survey of more than 16,000 
nearby F and G dwarf stars from the Nordstrom et al. (2004). Nevertheless, it is important 
to underline that all the observational data used in this study were taken from Holmberg et 
al. (2007), which brings an improved version of the previous catalogue, with new calibrations 
for all the relevant stellar parameters. All the selected objects are low-mass stars in the b — y 
range (0.2 — 0.6) corresponding to an effective temperature ranging from 4500 K to 7800 
K and a mass range 0.65Mq— ~ 2.0Mq. The low-mass limit at 0.65Mq reflects the red 
colour cutoff. Individual error estimation for all masses in the catalogue is about 0.05Mq. 
For a complete discussion on the observational procedure, calibration and error analysis the 
reader is referred to Nordstrom et al. (2004). 

The catalogue comprises a sample of 16,682 stars containing visual double stars and 
spectroscopic binaries of all kinds. The total number of binary stars of any type is 5,622. 
The sample thus comprises 11,060 stars that are not known to be double; of these, 7,817 
have measured radial velocities consistent with their being true single stars. 

The projected rotational velocities V sini are only given to the nearest km/s, and from 
30 km/s and upwards only to the nearest 5 or 10 km/s. Sense the great majority of stars 
have rotations below 20 km/s, for the present analysis we have thus selected 11,818 stars 
with Vsini < 30 km/s which we shall call the "original sample". 

The completness limit of the sample is a weak function of b — y. The distribution of 
the photometric V magnitude for all stars shows that the original sample begins to depart 
from completeness near V = 7.8. This is the value we adopt in this work. A total of 4,473 
stars have magnitude V < 7.8, satisfy the completness criterium, and we define these as our 
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The stellar sample 
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"complete sample". 

4. Results 

Both Tsallis and Kaniadakis generalized maxwellian are two parameter (q — o q and 
k — a K ) nonlinear functions given by (151). Since our aim is to test the deviation of the 
velocity distribution from the standard maxwellian we have fitted a g-maxwellian (Tsallis) 
and a K-maxwellian (Kaniadakis) distribution to the observed distribution using a nonlinear 
least squares fitting. We apply a nonlinear regression method based on the Levenberg- 
Marquardt algorithm. The value of the distribution parameters together with the reduced 
X 2 are given in Table 1. 

We can clearly see that the V sin % distribution of main sequence field stars in this sample 
does not obey a standard Maxwellian function since the values of q and K are significantly 
different from 1 and 0, respectively. 

Interestingly, the values of q and k do not change whether we consider the original 
sample or the complete sample. Nevertheless, the widths of the distributions are notably 
different. For the original sample a q = 4.25 ± 0.21 and a K = 4.95 ± 0.15, while for the 
complete sample we have a q = 5.08 ± 0.25 and o K = 5.92 ± 0.19. Even though the difference 
is significant, it is within the precision of the projected rotational velocities measurement, 
that is, 1 km/s. 

In the upper panel of Fig. [T] we show the best fits for the histogram of the observed 
distribution of V sin i according the results in Table 1 for the original sample. The q and k- 
Maxwellian functions are represented, respectively by the thin and thick lines. The dashed 
line represents the standard Maxwellian function. The distribution of observed V smi is 
without a doubt more adequately fitted by either a q and K-Maxwellian function. This can 
be more clearly seen in the lower panel of Fig. [T] where we have plotted the logarithm of 
the distribution divided by V sinz, that is log(<fr(y)/y), as a function of (V^sini) 2 so that the 
standard Maxwellian is represented by a straight (dashed) line. The Kaniadakis distribution 
is represented by the thick line while Tsallis distribution by the thin line. We observe that 
both non-gaussian distributions fit well the observed data, although the Kaniadakis function 
fits the data slightly better than Tsallis function. This can be seen in Table 1 where the 
values of the reduced-^ 2 is always smaller for the Kaniadakis fit compared with the Tsallis 
one. Also, the uncertainties in the parameters of the Kaniadakis distribution are smaller 
than in the case of Tsallis distribution. 

In Fig. [2] we present the fitting to the complete sample using the parameter given in 
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Table 1. The results are similar to the original sample. The q and ^-distribution give the 
best fit with Kaniadakis function being slightly better. 

Finally, we have found no noteworthy difference on the distribution of the projected 
rotational velocity between single and binary stars (see Table 1). Despite the fact that the 
values of q and a q and k and a K are slightly larger for binary stars, the difference are, in 
most case, within the statistical fluctuation. 

5. Conclusions 

In this work we have used non-gaussian statistics to investigate the observed distribution 
of projected rotational velocity of a magnitude complete sample of more than 16,000 nearby 
F and G dwarf stars. We have fitted a standard maxwellian and the q and /t-maxwellian 
based on the generalization of the power-law statistics proposed by Tsallis and Kaniadakis. 
We conclude that the V sin % distribution deviates significantly from a standard maxwellian 
and the best fit is attained for Tsallis distribution with q = 1.521 ± 0.043 and Kaniadakis 
distribution with k = 0.667 ± 0.033. These values correspond to all stars (single + binary) 
in the complete sample. 

As far as we are aware, this is the first time that Tsallis and Kaniadakis statistics are 
tested for such an exceptionally large sample of stars. This gives us an excellent degree 
of confidence on the fact that the distribution of stellar rotational velocity does not follow 
a Maxwellian-Boltzmann law as suggested by the pioneer studies by Chandrasekhar and 
Munch (1950) and Deutsch (1970). The present results show clearly that for the stellar rota- 
tion the non-extensivity holds, and that the distribution of the observed rotational velocity 
is explained by a generalized power-law statistics in the spirit of Tsallis and Kaniadakis 
statistical mechanics. 
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Table 1: Best values of the parameters of Tsallis (q and a q ) and Kaniadakis (k and a K ) 
distribution determined using \ 2 test for the rotational velocity of stars in the Nordstrom et 
al. catalogue. 



N 


9 


a-q (km/a) 


rcduced-x^ 




Ok (km/s) 


rcduccd-x^ 


Original Sample 
















all stars 


11818 


1.525 ± 0.041 


4.25 ± 0.21 


4493 


0.667 ± 0.029 


4.95 ± 0.15 


3088 


single 


6888 


1.484 ± 0.053 


4.13 ± 0.25 


422 


0.636 ± 0.038 


4.74 ± 0.17 


400 


binary 


4930 


1.551 ± 0.037 


4.63 ± 0.21 


223 


0.687 ± 0.031 


5.47 ± 0.17 


207 


Complete Sample 
















all stars 


4473 


1.521 ± 0.043 


5.08 ± 0.25 


550 


0.667 ± 0.033 


5.92 ± 0.19 


424 


single 


2382 


1.507 ± 0.046 


4.92 ± 0.26 


196 


0.655 ± 0.035 


5.69 ± 0.19 


142 


binary 


2091 


1.529 ± 0.049 


5.34 ± 0.31 


149 


0.672 ± 0.041 


6.26 ± 0.25 


133 



- 9- 




Fig. 1. — Observed distribution of the rotational velocity for all F and G field stars in the 
original sample. In the upper panel we show the histogram (number of stars N in each 
1 km/s bin as a function of Vsini), whereas in the lower panel, we show the logarithm of 
N/Vsmi as a function of (V sini) 2 . The curves represent the best fitted maxwellian (dashed 
line), Tsallis (thin line) and Kaniadakis (thick line) distribution. The fitting parameters are 
in Table 1. 
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Fig. 2. — Observed distribution (histogram) of the rotational velocity for all F and G field 
stars in the complete sample. In the upper panel we show the histogram (number of stars N in 
each 1 km / s bin as a function of V sin i) , whereas in the lower panel, we show the logarithm of 
N/Vsmi as a function of (Vsini) 2 . The curves represent the best fitted maxwellian (dashed 
line), Tsallis (thin line) and Kaniadakis (thick line) distribution. The fitting parameters are 
in Table 1. 



